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INTRODUCTION 

Study of stepped erystalline surfaees offers an exquisite 
combination of questions of fundamental physical inter- 
est and importance for technological development. These 
surfaces play a central role in the modern electronics in- 
dustry, where they serve as templates in the computer 
chip manufacturing process. The ability to model and 
quantify the evolution of such surfaces and ultimately 
engineer surface structures over a wide range of length 
and time scales will be essential for designing the next 
generation of computer chip components Even on 
the flattest of surfaces, e.g. 81(111), steps are inevitable 
on the submicron scale due to minor misorientations from 
facet directions. For metallic surfaces, steps appear on 
the scale of lOOnm. In some cases they provide a mech- 
anism of stress relief Q. 

In applications to electronics, template surfaces are of- 
ten grown epitaxially: material, usually in the form of 
atoms or molecules, is sputtered onto a crystalline sub- 
strate; the goal is to grow layer by layer with as few de- 
fects as possible. Many excellent books discuss this gen- 
eral problem from a variety of perspectives (e.g., [3, 0,13]). 
If one starts with a flat surface, growth begins at ran- 
dom nucleation sites. Islands expand around these sites. 
Eventually the boundaries of the islands meet. If the 
atoms in two abutting islands are in different domains, 
then a domain wall forms between them. Such walls per- 
sist during growth and cannot readily be annealed away. 
Instead, it is advisable to begin with a substrate that 
is intentionally slightly misoriented (typically by a few 
degrees) from flat, high-symmetry orientation. Such sur- 
faces are called "vicinal" since they are in the vicinity of 
the facet orientation. They consist of a series of atomi- 
cally flat terraces, separated from one another by surface 
steps — boundaries where the surface height changes by 
an atomic unit. Ordinarily the adsorption energy of a 
deposited atom is greatest at the crease at the lower edge 
of a step, since there it can bind to the largest number of 



other atoms. Thus, if the temperature is high enough so 
that the atoms diffuse relatively rapidly and the flux low 
enough, the atoms will be more likely to attach to the 
step edge than to meet another deposited atom to form 
a nucleation center on a terrace. Then, in this regime of 
"step-flow growth," the steps will gradually move across 
the terrace until, after deposition of a monolayer, the sur- 
face looks very similar to the initial vicinal surface, only 
with one more layer. 

Since steps play such a fundamental role, it is crucial to 
understand their properties and especially to clarify the 
few basic parameters that determine their behavior. The 
lowest-energy excitation of a stepped surface is the kink, 
a unit deviation perpendicular to the mean direction of 
the step. Since the energy to create an isolated atom or 
vacancy defect on a terrace is several times that to create 
a kink, kinks are the predominant defect on equilibrated 
surfaces at low temperatures. While these kinks cost en- 
ergy, they contribute to the entropy in the usual way, so 
that in equilibrium this competition leads to the step free 
energy per length /?, or line tension [lo| . In an expansion 
of the projected free energy per area of a vicinal surface 
(i.e.; the surface free energy per area of the vicinal surface 
projected onto the terrace plane), (3 is the coefficient of 
the density of steps. (At the roughening transition of the 
terrace plane, /3 vanishes and the projected free energy 
per area becomes quadratic in the density of steps; the 
vicinal surface is in this technical sense rough already at 
the lower temperatures under consideration.) In other 
words, the line tension indicates the extra energy associ- 
ated with a unit length of a step. 

Steps need not run along high-symmetry directions, i.e. 
the surface normal (or azimuthal misorientation) can be 
at an arbitrary polar angle 9. Furthermore, the border 
of a single-layer island (or vacancy island) is just a step 
that is a closed curve — somewhere between circular and 
polygonal — rather than a nearly straight line. The cele- 
brated Wulff construction [l]| uses f3{6) to determine the 
equilibrium crystal shape that minimizes the free energy 
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at constant area. 

The thermal excitation of kinks along the step leads 
to meandering. Such fluctuations are constrained by the 
stiffness (3(9) = (3{9) + f3"{9), which weights the squared 
slope of the step relative to its mean direction 0, H, Q . 
(There is an analogy between the ensemble of spatial con- 
figurations of the steps in 2D (two dimensions) and the 
world lines of particles evolving in ID; in this picture, 
the deviation of the slope is analogous to velocity, so the 
wandering corresponds to kinetic energy, with stiffness 
playing the role of mass.) Furthermore, it is the stiffness 
that weights the curvature in the Gibbs-Thomson term, 
the curvature contribution to the step or interface energy. 
Thus, in many ways (3{9) is more fundamental than P{9) 
and in some situations is better defined [l^ . 

Accordingly, stiffness is one of the three parameters of 
the step continuum model, which retains distinct steps 
but coarse-grains them into continuous "strings." The 
other two parameters are the strength of the inverse- 
separation-squared repulsion between steps and the char- 
acteristic rate of the mechanism dominating step kinet- 
ics. Then, at the mesoscopic scale, such surfaces can 
be envisioned as a collection of steps separating atomi- 
cally fiat terraces and effectively tracing out surface con- 
tour plots. By tracking the net movement, fluctuations, 
and interactions of steps, the evolution of the entire sur- 
face can be monitored. The actual motion of atoms — at 
much smaller length and time scales — which underlies 
this behavior appears only in the parameters, which can 
be estimated from either calculation (or measurement) of 
atomic processes or by fitting to data at the mesoscale. 
This idea, at the heart of the step-continuum model Q , 
turns out to be an extremely efficient way of tracking 
surface evolution over a wide range of experimentally 
relevant length and time scales. This picture has been 
used to successfully account for a broad range of step 
properties, such as mound decay, step and island fiuctu- 
ations, cluster diffusion, electromigration, and ripening 
. Since one can never be certain of including all signif- 
icant atomistic-scale processes, the most meaningful test 
of the step continuum model is the check that the same 
set of parameters describes quantitatively all these var- 
ied phenomena. One can then return to the microscopic 
picture, invoke a simple model, and fit the key model pa- 
rameters as effective values to reproduce the mesoscopic 
behavior. If reliable energetic calculations give agreement 
with these predictions, all the better. 

In this paper we concentrate on the two densest, 
highest-symmetry faces of an fee crystal, namely {100} 
and {111}, which have square and hexagonal symmetry, 
respectively. (For {111} surfaces, if one imagines adsorp- 
tion into both kinds of 3-fold sites [fee and hep] , one has 
a honeycomb. This feature is unimportant for what we 
consider here.) Late transition and noble metals have fee 
crystal structure and make good substrates for the sorts 
of experiments envisioned here, since they are relatively 



soft, with atomic motion occurring adequately for equili- 
bration at room temperature or at somewhat higher. As 
the close-packed Bravais structure, there is less angular 
dependence on bonding, making near-neighbor models 
better approximations than for bee metals. Note that 
{100} and {111} faces of all cubic Bravais crystals have 
square and hexagonal symmetry, respectively. So does Si 
and other systems with the diamond structure. 

If one assumes that step adatoms interact with 
only nearest-neighbors (NN) or next-nearest-neighbors 
(NNN), then it is possible to derive exact solutions for 
the line tension based on the Ising or sohd-on-solid (SOS) 
models. These solutions are implicit, however, making 
their implementation into numerical simulations time- 
consuming and computationally demanding, particularly 
when dealing with the stiffness, which requires two ad- 
ditional derivatives of the implicit line tension. For sim- 
plicity, then, numerical studies often 1^, 31 (though by 
no means always [3l) assume an isotropic line tension 
and stiffness. Except at high temperatures where an is- 
land structure is nearly circular, this approximation turns 
out to be poor, especially near facet orientations, where 
the line tension is notably smaller, and would be much 
smaller if there is a quasi-straight-edge (2D "facet" — at 
macroscopic scales the island is rounded at finite tem- 
perature [3]). For the stiffness, the problem is more 
severe, since (3" [9) is large near these special directions 
(though not infinite as it would be in 3D), leading to 
a stiffness much greater than the value at general ori- 
entations. Since it is easier to compute (3{9) in such 
high-symmetry directions, such values have been used 
to characterize (3{9) at general 0, thereby considerably 
overestimating the typical stiffness. 

The next simplest approximation assumes a sinusoidal 
variation reflecting the substrate symmetry 17[. Again, 
there are shortcomings to this procedure, especially near 
facet orientations, about which polar plots of the line ten- 
sion as a function of angle reveal sharp, cusp-like minima, 
implying the stiffness blows up (since the curvature of a 
cusp is infinite). Such temperature- independent simpli- 
fications preclude quantitative comparisons with experi- 
ment [13]. 

In this paper we construct expressions for [3{9) and 
I3{9) that are well behaved analytically, being continu- 
ous and twice differentiable and that give an accurate 
accounting at all orientations and relevant temperatures. 
For convenience, all derived formulas are summarized in 
Table 5.1. While not especially simple, they are straight- 
forward to construct and easy to implement in numer- 
ical codes such as used in finite-element investigations 
Til . 20| . making quantitative comparisons with dynamic 
experiments possible. We thus expect our results to be 
widely applicable. 

Our approach begins with simple, low-temperature for- 
mulas for the orientation dependence, on face-centered- 
cubic (fee) surfaces, of the {001} and {111} stiffness and 
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line tension that we derived in two recent papers 2^, 2^ . 
(This approach is rooted in the lattice-gas perspective, 
so is complementary to Shenoy and Ciobanu's study of 
stiffness anisotropy based on elasticity theory [l^]-) Our 
formulas assume the step fluctuations are dominated by 
the rearrangement of geometrically forced kinks — kinks 
that are not thermally activated. At temperatures low 
compared to the surface roughening temperature (for 
noble metal surfaces, such as Ag and Cu, room tem- 
perature is considered "low"), the formulas only fail for 
steps having a negligible number of forced kinks; that is, 
steps oriented very close to the high-symmetry direction. 
When the step angle is exactly 0° (aligned with the high- 
symmetry direction), the formulas predict a cusp in the 
line-tension and an infinite step stiffness. Here we cor- 
rect for the non-analytic behavior by splicing our simple, 
low-temperature formulas with small-angle expansions of 
the exact, implicit solutions based on the Ising and SOS 
models. 

In the following section, we describe the details of a 
general expansion for the stiffness and line tension that 
is continuous and twice-differentiable. In sections III and 
IV, we apply this expansion to fee {111} and {001} sur- 
faces, respectively, to derive surface-specific formulas for 
the stiffness and line tension. In the final section, we offer 
concluding remarks as well as a synopsis of the derived 
expressions. 

EXPLICIT ANALYTIC APPROXIMATION 

At the microscopic level, the step stiffness and line ten- 
sion arise from the energy and rearrangement of step edge 
kinks. It is therefore natural to decompose (3{d) and f]{d) 
into two contributions: one part originating from geomet- 
rically forced kinks and one part from thermally activated 
kinks. Geometrically forced kinks, depicted in the inset 
of Fig. [H are present at all temperatures, and give the 
step an overall orientation 9. The further 6 is from the 
high-symmetry direction, the greater the number of ge- 
ometrically forced kinks. Thus, at lower temperatures, 
as long as the orientation angle of a step is greater than 
some small, temperature-dependent cross-over angle 9c, 
there are many geometrically forced kinks and relatively 
few thermally activated kinks, suggesting f3{9) and (3{9) 
can be well described by formulas based on geometrically 
forced kinks alone. 

As an example, we have recently derived a remark- 
ably simple, low-temperature formula for the {111} step 
stiffness assuming only NN adatom interactions and ge- 
ometrically forced kinks: 



fcsT _ sin(36i) 
W) ^ 2V3 ■ 



(1) 



angles, but predicts an infinite stiffness when 9 = 0. For- 
tunately, the exact, implicit solution based on the NN 
Ising model can be explicitly written for steps having this 
orientation. We can therefore expand the exact solution 
about 9 = and splice it with our low-temperature so- 
lution at 9c, thereby producing an explicit form for f3{9) 
valid at all angles. This idea is illustrated in Fig. ([1]). 
Here, an additional orientation-independent contribution 
to the stiffness from thermally activated kinks A is also 
included for completeness. Similar to high-symmetry 
steps, the stiffness of maximally kinked steps {9 = tt/6) 
can be exactly obtained from the NN Ising model, so that 
A can be determined explicitly. 

To generalize this approach, we assume f3{9) and (3{9) 
are well described at angles greater than 9c by simple, an- 
alytic functions representing contributions from geomet- 
rically forced kinks. Explicit forms for these functions 
[2^ . [2^ will be discussed later. For now, to be general, 
we simply write them as f{9). 

At sufficiently low temperatures, 9c is small, so we may 
accurately represent l3{9) and the inverse stiffness P~^{9) 
at angles less than 9c using small- angle expansions. (We 
expand the inverse stiffness because, in the 9=0 limit, it 
vanishes at low temperatures, making it mathematically 
better behaved than the stiffness itself, which diverges). 
Specifically, we construct an approximant X{9) to rep- 
resent the dimensionless form of the function we wish to 
expand — either (3{9)a\\/ [ksT) or ksT / {I3{9)a\\), where 
a|| is the close-packed distance between atoms (i.e. the 
atomic diameter), and /csT is the Boltzmann energy — we 
define 



X{9) 
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ii9 <9c 
ii9> 9c 
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To fully specify this function, we must find the appro- 
priate expansion coefficients, a„. We obtain their values 
by matching Eq. ([2]) and its higher order derivatives with 
the exact solutions at 9=0 (which can be systematically 
obtained) and the approximate (yet accurate) solutions 
obtained from f{9) a,t 9 = 9c, analogous to performing a 
spline fit [2^. Specifically, for the boundary conditions 
at 9 = 0, we have 

„ _ d^x{o) 



n< N 



(3) 



where dJ}X{0) = d''X{9)/d9''\e=o. The remaining N 
coefficients are found from the boundary conditions at 
9 = 9c, which form a set of N coupled linear equations: 



2N-1 
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(4) 



At sufficiently low (but experimentally relevant) tem- 
peratures, the formula works well for steps at nearly all 



where m = 0, 1,2, ...,N - 1. 

For use in continuum models, (3(9) should be continu- 
ous and twice-differentiable. To ensure the second deriva- 
tive remains continuous at 9 = 9c, this requires, at min- 
imum, A'^ = 3. In this case, Eqs. ^ are simultaneously 
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FIG. 1: (Color online) The contributions to the step stiffness can be decomposed into parts originating from geometrically 
forced kinks (lower blue region bounded from above by the line labeled "low-T" ) and thermally activated kinks (the remaining 
red region, bounded from above by the line labeled "exact"). At relatively low-temperatures, the {111} step stiffness is well 
approximated at angles greater than Qc by a relatively simple, explicit function /(6'), since the thermal part is evidently 
insensitive to angle. To account for all angles, the formula can be spliced with a small-angle expansion of the exact NN Ising 
model solution (from which explicit forms for the stiffness can be obtained at S = and at 7r/6, depicted here by hollow circles). 
The solution at 7r/6 is used to determine A. The expansion coefficients a„ are obtained by matching the solutions at S = 
and Qc- The inset depicts a step edge from above. Each square represents an adatom which is part of the step edge. The 
upper-most square represents a thermally excited adatom, which forms four thermally-activated kinks. The remaining kinks 
are geometrically forced — they must be present to give the step edge an overall angle Q. 



solved to give: 
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as 



20(/ - X) ~ 8/' + (/" - 3X") Bl 
2 613 

-30(/ - X) + 14/' - {2f' - iX 

i2(/ -X)- 6/' + if" - X") 



(5) 



(6) 



(7) 



where the prime represents differentiation with respect to 
6; for brevity we write / = f{6c) and X = X{Q). Note we 
have also used Eq. which implies oo = X, ai = X', 
and a2 = X" /2. Because both the line tension and the 
stiffness are continuous and symmetric about 6=0, we 
know that ai = X' ~ 0. In the remaining sections we 
apply this approximation to specific cases where explicit 
forms for X and / can be obtained. 



{111} SURFACES WITH NN INTERACTIONS 

For {111} surfaces with only NN adatom interactions, 
Zia found an implicit form for the full orientation depen- 
dence of the step line tension [2^ : 



knT 



MS)iJiiO,T/T,)+fj^{d)M0,T/T, 



(8) 



where 770(6') = (2/^3) sin(6'), 7]±{e) = cos{e) ± 
(l/\/3) sin(6l). Here Tc is the critical temperature of 



the NN lattice-gas model. The ■0's are solutions of the 
pair of simultaneous equations for the angular constraint. 



sinh(i/)i - i-02) cosh(i?/'2) 
sinh(i/)2 — ^"01 ) cosh(i?/;i) 

and the thermal constraint, 

cosh')/'! -I- coshi/)2 + cosh(V'i — 1P2) 



r]o_ 

V- 



(9) 



(10) 



where y = y/{3z + l)/z{l - z) and z = 3"^-/^. The 
latter can be rewritten z = exp(— 2efe/fcBT), where is 
the energy of a kink on a close-packed step and 



From Eqs. ([OTTOl) it follows that 



(11) 



0i(O) = ^V'2(O)-cosh-i(^^). (12) 

With ■01(0) and V'2(0) in hand, wc can differentiate the 
constraints, Eqs. (|9ll0p . set = 0, and systematically 
solve for all the higher order derivatives of the -i/i's, which, 
according to Eq. (I8|), are sufficient to find the higher or- 
der derivatives of (3. We will utilize these higher order 
derivatives to derive explicit, analytic approximations for 
the stiffness and line tension. 
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FIG. 2: (Color online) In the left plot, the orientation depen- 
dence of our explicit approximation for the {111} step stiffness 
(solid lines) and its inverse (inset, solid lines) are compared 
to the exact, implicit solutions (symbols). Because of the six- 
fold symmetry of the solution, only the positive half of the 
first sextant is shown (the negative half is mirror-symmetric) . 
The right plot shows the values used for dc (solid dots) in 
the construction of the left figure and the corresponding ex- 
ponential fit (solid line), good over the temperature range of 
interest. The fit is expressed in terms of the kink energy et, 
which is related to Tc by Eq. (|ll|l . The inset shows the sum 
of squared vertical deviations (^ x^) versus angle in a least 
square fit for 9c. At each temperature, 6c is the angle that 
minimizes this sum. 



the first term has any 6 dependence, /' and /" are simple 
to calculate. 

Now only X and its first two derivatives need to be 
determined. As mentioned in the preceding section, these 
can be systematically determined. In particular, we find 
|see Eq. (23) for a derivation of X in our earlier paper 



X 
X' 
X" 



3(y - 1) 



= 0, 



y3 - 2?/2 - I5y + 36 



2(j/-i)v/y'-2y-3 



(14) 
(15) 
(16) 



Of course, based on symmetry, we already knew that 
X' = 0. 

By combining the functional forms for / and X and 
their derivatives with Eqs. ^E^, we can plot the stiff- 
ness and compare it to the numerically evaluated exact 
solution. We show this comparison in Fig. [2l where 6c 
was determined at a variety of temperatures by doing 
least square fits to the exact solution. The agreement 
shown in Fig. [2] is very good at low-temperatures and is 
quite reasonable at temperatures all the way up to Tc/5. 
(This behavior is remarkable since slightly above Tc/5. 5, 
dc becomes greater than 30°, i.e., the power series is used 
for the entire range of orientations. Once \6c\ > 30°, the 
slope of ksT / {[3{9)a\\) no longer vanishes at 30°.) At 
higher temperatures, the angular dependence becomes 
negligible, so [3{6) become isotropic. 

The right plot in Fig. [2] shows the values used for 9c, 
along with an exponential fit: 



Step Stiffness 

In this case, X{6) = ksT / {(3{9)a\\), which is six-fold 
symmetric for {111} surfaces with only NN adatom inter- 
actions. To utilize our explicit analytic approximation, 
we require f{9) — the contribution to the reduced stiffness 
from geometrically forced kinks — which, in the first sex- 
tant (— 7r/6 to 7r/6), takes a relatively simple form [23[: 

f{9) = { sin(30) + ^^y"" , _ 1 1 . (13) 
2^/i\ - 10y2 + 9 ; 

The last two terms, called A in Fig. [TJ are included to 
ensure f{9) matches the exact solution for steps with 
orientation angle 9 = tt/Q. The physical origin of the A 
terms is the thermal fluctuations of a maximally kinked 
step. Such fluctuations are relatively inexpensive in 
terms of energy. They dominate the fluctuation contri- 
bution while a significant fraction of the step is not close- 
packed, so that the thermal contribution for such orienta- 
tions is relatively independent of orientation. Since only 



9c{T) « 11.2 eM-^k/ksT) = M2[°]{^)-'^^-'^ . (17) 

The second form uses the units in Fig. [21 reexpressing 
the prefactor in degrees and the exponent in Tc/T . The 
Arrhenius decay reflects the importance of thermally- 
activated kinks for \9\ < 9c- 

Step Line Tension 

We follow the same procedure for the line tension. In 
this case X[9) = (3[9)a\\/kBT . Corresponding to the 
T = divergence of the stiffness at 6* = is a cusp in the 
line tension, indicating a facet in the equilibrium shape. 
At finite T the cusp, like the divergence, vanishing since a 
facet on a 2D structure corresponds to ID long-range or- 
der. Specifically, the contribution (in the first sextant) to 
the line tension from geometrically forced kinks is fairly 
simple [iil: 

f{9) — — ?/+ In z — r^-i- In r/-|_ + In r/_ + ?/o In ?/o- (18) 
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FIG. 3: (Color online) The orientation dependence of the ex- 
plicit approximation for the {111} line tension (solid lines) is 
compared with the numerically evaluated exact result (sym- 
bols). Because of the six- fold symmetry, only the positive half 
of the first sextant is shown. (The negative half is mirror- 
symmetric.) 



Just as for the stiffness, we systematically determine X 
and its first two derivatives by differentiating the exact 
solution, Eqs. PfH)) . 



X 
X' 
X" 



^=2cosh-^^-l 



2y v/y2 - 2y - 3 
3(y-l) 



-X. 



(19) 
(20) 
(21) 



The last equation can be rearranged to find the reduced 
stiffness at 6* = 0, as expressed earlier in Eq. With 
these parameters in hand, we compare our approximation 
for the full orientation dependence of the reduced line 
tension with the exact, numerically evaluated solution in 
Fig. m For the critical angle, we use Eq. (jl7|l . As before, 
the fit works remarkably well at temperatures as high as 
Tc/5. 



{001} SURFACES WITH NN AND NNN 
INTERACTIONS 

For {001} surfaces with just NN interactions, an ex- 
act, explicit form for the full orientation dependence of 
the line tension was first determined by Abraham and 
Reed [13]. For such surfaces, however, NNN interactions 
are often significant [l^l, so it is desirable to find a so- 
lution including their effects. We denote by R the ratio 
of NNN to NN adatom interaction strengths; the latter 
is assumed to be attractive (negative), so a positive R 
indicates that the NNN interaction also is. 

Although no exact solution to the Ising model with 
both NN and NNN interactions exists, the solid-on-solid 
(SOS) model provides an excellent approximation at rea- 



sonable temperatures (^ Tc/2 based on our comparisons 
with the imaginary path weight random-walk method de- 
veloped by the Akutsus [1^ ) . This model can be solved 
exactly |22l |. yielding the following implicit form for the 
reduced line-tension: 



knT 



[p{e) sme + g{p{e))]cost 



(22) 



where p{6) is found by inverting 

2 sinh p sinh S 



tan^ 



(cosh 5* — coshp) [2 sinh — (coshS* — coshp)(y + 1)] ' 

(23) 

while g{p) is 



g{p) = - In 



y 



1 



sinh S 



y — 1 1 — y cosh S — cosh p 



(24) 



Here y = I - 2z^, S = -(i? + 1/2) Inz, z = (1 
^^-2Te/T ^ exp(-2efc/fci3r), while is the critical 
temperature for i? = (just NN interactions): 



ksTc 



ln(l + V2), 



(25) 



where the kink energy now refers to a close-packed 
step on an {001} surface. Wc will utilize the exact, im- 
plicit solution Eqs. (|22ll24p to determine the parameters 
required to find an explicit approximation for the stiffness 
and line tension below. 



Step Stiffness 

To begin, we let X{9) = keT / {j3{e)a\\). The symme- 
try of {001} surfaces require X{9) be four-fold symmet- 
ric. Accounting for just geometrically forced kinks, the 
reduced inverse stiffness is well approximated in the first 
quadrant (—7r/4 to 7r/4) for \9\ > Oc by the following 
function 12211: 



/W = ^v/l-ysm(2^). 



(26) 



By differentiating Eq. ([26|) , /' and /" are easily obtained. 

To determine A, A', and X" (and, potentially, any 
higher order derivatives), we utilize the exact solution of 
the NNN SOS model. Eq. (QSl), for example, imphes that 
po = when = 0. With some effort, it can be shown 
that 

X = 2 sinh g 

(coshS'-l) [2sinh5-(cosh5'-l)(y-Hl)] ^ ' 
A' = (28) 



A" = 



1 2 cosh 5*+ 1 
A coshS'-l 



coshS'— 1 y+l 



sinh 5* 



A 



.(29) 



As required by symmetry. A' = 0. 
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FIG. 4: (Color online) In the upper-left plot the orientation dependence of the explicit approximation for the {001} step stiffness 
(solid lines) and its inverse (inset, solid lines) are compared to the exact, implicit solutions (symbols). Because of the four- fold 
symmetry of the solution, only the positive half of the first quadrant is shown (the negative half is mirror-symmetric). The 
upper-right plot shows the values used for Oc (solid dots) in the construction of the upper-left figure and the corresponding 
exponential decay fit (solid line) good over the temperature range of interest. The fit is expressed in terms of the kink energy 
Efe which is related to Tc by Eq. (|25p . The inset shows the sum of squared vertical deviations (X]x^) versus angle in a least 
square fit for Oc- At each temperature, 9c is the angle that minimizes this sum. The two lower plots show the {001} inverse 
stiffness for a variety of different R at two temperatures, Tc/9 and Tc/5 (the extremum of the temperature range of interest). 
Notice that for a given temperature, all curves align at an angle greater than the largest critical angle ff^""^ . This behavior 
means Qc, practically speaking, does not depend on R at these temperatures. 



Combining the functional forms for /, X, and their 
derivatives with Eqs. ([2][7l), we can plot the inverse step 
stiffness and compare it to the numerically evaluated ex- 
act solution, just as before. We show this comparison in 
Fig. m where was determined by doing least square 
fits to the numerically evaluated exact solution (with 
R = 1/5). The agreement shown in Fig. [4] is excellent 
at low-temperatures and is very reasonable at tempera- 
tures all the way up to Tc/5, as was the case for the {111} 
solution. 

Although it was not initially obvious, the relative size 
of the NNN interaction R has little effect on 9c. This 
fortuitously implies that a single 9c works for all values 
of R, as depicted in the lower plots of Fig. [D 

With this in mind, the values used for dc were deter- 
mined just as they were for the {111} case, but with 
R = 1/5. These are shown in the upper-right plot of 
Fig. m as well as a simple fit that is accurate over the 
temperature range of interest: 



9c(r) « 6.72 e-<^'=/'=«'^ = 385[°](1 + V2y 



(30) 



As for Eq. (fTT)) , the second form uses the units in Fig. |4l 
reexpressing the prefactor in degrees and the exponent 
in Tc/T . Again, the Arrhenius decay is anticipated since 
9c represents the angle below which thermally activated 
kinks on close-packed segments become important. 

Finally, we point out that the {001} step stiffness is 
much more anisotropic than its {111} counterpart. In 
fact, at Tc/6 the anisotropy is as large as the {111} 
anisotropy at Tc/Q. Furthermore, 9c is less sensitive to 
temperature than its {111} counterpart. This follows 
from the relative ease of thermally activating kinks on 
{111} steps, requiring only the breaking of one NN bond, 
as compared to two for {001} steps. For {111} steps, 
then, the angle 9c below which thermally activated kinks 
become important is larger than for {001} steps. 

However, there is no need to include an "off-angle" cor- 
rection A as was needed for the {111} case, at least in the 
case of just NN interactions (i?=0, y = — 1). In that case 
one can readily find the difference between the stiffness 
of the exact result and /(7r/4) from Eq. (US]): 
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1 

^ L\/l-4scch2(e;,/fcsT)tanh2(efc/fcBT) 

2V2sinh"(£fc/fcBr) 
cosh(2efc/fci3r) -3' 



(31) 



Over the range of temperatures of interest here, numeri- 
cal evaluation shows Aioo is negligible. 



Step Line Tension 




We proceed as usual, letting X{9) = (3{9)a\\/kBT . The 
contribution from geometrically forced kinks is found by 
solving the low-temperature form of Eq. (|23p . which be- 
comes quadratic in ef^^. Solving gives 



— ysai[29) + y sin 6 — cos ( 
(1 +y)sine 



(32) 



Plugging this into Eq. (f22|) yields an excellent approxi- 
mation f{6) for the reduced line tension X{9) valid in 
the first quadrant (-7r/4 to 7r/4) for \9\ > 9c: 



f{9) = cose* 



+ In ■ 



(1-?;) sin 6i-hcos 61-^1- ysin(26') 



(1 + y) (sin6'-cos6'+^l-ysin(20) 



- sin 6* 



S' + ln 



^yT^^^y^n(29) + y sin 9 — cos ( 
(1 + y)sm9 



FIG. 5: (Color online) The orientation dependence of the ex- 
plicit approximation for the {001} line tension (solid lines) is 
compared with the numerically evaluated exact result (sym- 
bols). Because of the four-fold symmetry, only the positive 
half of the first quadrant is shown (the negative half is mirror- 
symmetric) . 



line tension and compare it to the numerically evaluated 
exact solution. We show this comparison in Fig.[5l where 
9c was determined from Eq. ([30|) and R = 1/5 (other 
values yield equally good agreement). As before, the ap- 
proximation works well at temperatures up to Tc/5 (and, 
in this case, perhaps even higher). 



Differentiating twice straightforwardly gives /' and /". 
Eq. (|33p can be written more compactly by defining 

and inserting w{9,y) = cos 9 — -^/l — ysm{29) /sin0, 
as done in Table 5.1. 

This leaves X and its derivatives. They too can be 
explicitly determined from the exact solution. Setting 
both 6* = and p (as Eq. ^ demands) in Eq. (|22ll . 
we find X: 



X 



.9(0) 
5 -In 



y 



1 



sinh S 



2/ — 1 1 — y cosh 5—1 
Similarly, it can be shown that 
X' = 0, 



(34) 



(35) 



X" = 



(coshS'-l)[2sinhS'-(cosh5'-l)(?/ + l) 
2 sinh 5 



-X. (36) 



This last equation can be rearranged to give the reduced 
step stiffness, as previously written in Eq. ([27|) . 

By combining the functional forms for / and X and 
their derivatives with Eqs. ^i7^, we can plot the reduced 



(33) 



SUMMARY AND CONCLUDING REMARKS 

We have constructed explicit, twice-differentiable ap- 
proximants for the full anisotropy of step stiffness and 
line tension on both {001} and {111} surfaces of fee crys- 
tals, the metallic systems that have been subjected to the 
greatest scrutiny with regard to island properties [30l.l3lt. 
Our expressions are accurate over a broad range of ex- 
perimentally relevant temperatures; they fail only when 
the stiffness is nearly isotropic, i.e., when their use is no 
longer required. Implementation into continuum simu- 
lations is straightforward and efficient. They are much 
more usable than numerically extracting solutions from 
the underlying 6th-order equations, and more fiexible and 
convenient than constructing immense look-up tables as 
functions of angle and temperature from such a proce- 
dure. Our expressions are greatly superior to conven- 
tional explicit formulas for step stiffness and line ten- 
sion, which usually take the form of simple sinusoidal 
variation that neither carry temperature dependence nor 
accurately capture the anisotropy (extreme for the step 
stiffness) observed at lower temperatures. For clarity and 
convenience, we summarize our results in Table HI 

We have implemented these formulas into state-of-the- 
art finite-element simulations, where we have tested them 
by numerically determining the equilibrium shape of two- 
dimensional islands at various temperatures Sj] . Explic- 
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EXPLICIT APPROXIMATION FOR STIFFNESS & LINE TENSION 

( 2JV-1 - V _ 20(/-X)-8/' e,+(/"-3X") 

xie) := < 



/(e), e > 



ao 


= X 


as 


ai 


= 


a4 




_ X" 




02 


2 


as 



-30(/-X) + 14/' 9c-(2/"-3X") 



2 

12(/-A-)-6/' 6)c+(/"-A-") 6)g 

2~e; 



{111} Surfaces with NN Interactions 



Stiffness {X{e) « /cBT/an/?) 



Line Tension {X{e) ^ au^/k^T) 



X 
X" 

an 



3(i;-i) 



i;^-2y^-15i;+36 



2(!;-i)Va^-2!/-3 



- 1 



2cosh-i(^) 

3(H-1) 

—7?+ ln(2r;+) + r)^ Inr)- + r;o In 770 



X 

X" 



2 sinh 5 

(cosh S-l)12sinhS-(cosh S-l)(y + l)| 

J_ 2 cosh S + 1 _ A [ coshS-l y + 1 , y-] 
X coshS-1 L sinhS 2 J 



S - In 



y + l I 2 sinh S 
y~l 1 — y coshS — 1 



(cosh S-l)[2 sinh S- (cosh S-l)(y + l)] 
2 sinh S 



(l~y)(l + ^.(9,y)) 

™(e,y)) 



cose[5+lnii^ 
+ sine I^S + ln^i:^! 



-X 
+ 



y-w{0,y) I t 



{001} Surfaces with NN and NNN (= RxiViV) Interactions 

[°]=385 e-^B^fe^ ^^g-2/3Befe^(i + ^)-^^ S = (l + 2i?)/3Befe, y = l-2z« 



*ri± = cos 6* ± sin 6^, rjo = -7= sin 



\/3 ■ 



'j{6, y) = cote- CSC 6'^l-ysin(26l) 



TABLE I: (Color online) Summary of results for approximants of dimensionless inverse stiffness and line tension. X = X{Q), 
while / = f{Oc); /3b = (/cbT)^^, the subscript needed to distinguish it from the line tension. The upper part of the table 
(dark red) refers to the steps on the hexagonal-lattice face, with just NN interactions. The lower part (blue) refers to the 
square-lattice face; by setting i?=0, one retrieves the simpler formulas for just NN interactions. 



itly, this is done by finding the shape that minimizes the 
chemical potential (proportional to the product of the 
step curvature and stiffness) of a single, island-bounding 
step. We have also used these simulations to model is- 
land relaxation from an arbitrary shape to the equilib- 
rium shape and are currently using them to model re- 
cent experiments monitoring the relaxation of depinned 
Ag(lll) steps [32I. 
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